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I. INTRODUCTION 



Recent success in developing experimental methods for dripline nuclei, that in partic- 
ular allow exploration of halo phenomena in light nuclei, has put on the agenda a need 
for appropriate theoretical methods which take into account the peculiarities of weakly 
bound and spatially extended systems. For Borromean two-neutron halo nuclei ( 6 He, 11 Li, 
etc.) an understanding of the essential halo structure has been obtained in the framework 
of 3-body models [|IJ . Reactions involving these nuclei present however, at least a 4-body 
problem. The direct solution of 4-body systems is extremely difficult, and approximate 
methods are required. For high energy elastic scattering and relativistic fragmentation of 
Borromean halo nuclei, a 4-body Glauber method has been developed P,f| . For Coulomb 
breakup or electromagnetic dissociation (EMD) the first order (Alder- Winther) pertur- 
bation theory or an equivalent semiclassical treatment || has been used, but with exact 
3-body continuum wave functions ]5|-[7||. Also for the (anti)neutrino induced reactions on 
6 Li populating the 6 He and 6 Be 3-body continua, have proper final state wave functions 
recently been used ||. 

The most reliable information on properties of halo nuclei, especially for the low-lying 
part of excitation spectra, is experimentally obtainable by intermediate energy elastic and 
inelastic scattering and charge-exchange reactions. The distorted wave theory is the most 
common way to analyse such processes ||, but for halo systems their spatial granularity 
as well as peculiarities of their quantum structure have to be taken into account. The 
3-body interaction dynamics defines the low-lying part of excitation spectra, in particular 
the soft modes of Borromean systems, and has to be treated properly. Until now, only 



the hyperspherical harmonics (HH) method [10] is able to provide a formulation of the 
scattering theory to the 3-body Borromean continuum. The Faddeev equations technique 
|TT| has been developed to investigate breakup of 3-nucleon systems, but has hitherto not 
been applied to investigate the continuum in Borromean nuclei. The coordinate complex 
rotation method [0,0] gives Gamow and not scattering states, and is mostly suitable for 
searching for resonance positions or poles in the complex energy plane. 

In previous studies []],|n]], the HH method gave a very successful and comprehensive 
description of data on weak and electromagnetic characteristics of the 6 He and 6 Li sys- 
tems, and of the absolute values of differential cross sections of (p,n), (p,p') and (n,p) 
reactions to the bound states of A = 6 nuclei. These nuclei still represent the best test- 
bench for quantitative calculations for Borromean halo nuclei. In the present work we 
develop distorted wave theory for inelastic and charge-exchange reactions leading to the 
3-body continuum. For the continuum excitations of 6 He we perform a detailed analysis 
of inclusive excitation and differential cross-sections for beam energies in the range of 
the GANIL facility, where such experiments are in progress. Investigations of continuum 
spectra of 6 He are also subject of future experiments at Kurchatov Institute (Moscow), 
NSCL (Michigan), RIKEN (Tokyo), JINR (Dubna) and GSI (Darmstadt). 



II. SHORT PREAMBLE 

The known spectrum of 6 He contained until quite recently only the + bound state, 
the well known 2 + (1.8 MeV) 3-body resonance and then a desert in the 3-body a+n+n 
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continuum up to the 3 H + 3 H threshold at about 13 MeV JT^|. With radioactive nuclear 
beam techniques and dynamic approaches to 3-body continuum theory |10| new possibil- 
ities are opened, and we may now ask to what extent our knowledge of 6 He is complete, 
and what specific influence the halo has on the continuum structure. The so-called soft 
dipole mode suggested in [0,0 was the first example of this quest. That it is not a 
simple binary core-point dineutron resonance in n Li seems now widely accepted, nor is it 
in 6 He according to our recent calculations [Tj|. Is it rather a genuine 3-body resonance 
or just a dynamically induced very large dipole moment, or the consequence of two-body 
final state interactions? Soft modes of other multipolarity have also been theoretically 
suggested BIJ. Their presence needs clarification, both theoretical and experimental. 



The most natural way to observe soft modes in exotic nuclei is by inelastic excitation 
of radioactive beams. In the 6 He case, however, only results of fragmentation experiments 
without reconstruction of inclusive spectra have been published |p2|-f2~4jj. Other ways 
include transfer reactions (like 7 Li(n,d) 6 He at E n = 56.3 MeV ||25|| ) or charge-exchange 
reactions of (n,p) type on 6 Li. At E n = 60 MeV the 6 Li(n,p) 6 He reactions have been 
measured, but with poor statistics and limited angles ||26| . In heavy- ion charge-exchange 
reactions 6 Li( 7 Li, 7 Be) 6 He a broad bump at excitation energy ~ 6 MeV was observed 
2"7|j2"g|] , but different assignments were made about its multipolarity. 



In our recent work |T9|,20| we have used several methods to investigate the internal 
structure of the 3-body continuum, as well as the transition properties for accessible 
nuclear reactions in terms of nuclear and electromagnetic response functions. Our meth- 
ods have the advantage that, even off a resonance, the continuum structure can still be 
investigated while taking into account all final state interactions. 

We have predicted [|19| surprising richness of the 6 He continuum structure, and applied 
elaborate methods to explore the nature of different modes of excitation of the Borromean 
halo continuum ||20|| . In 3-body dynamics we have found two kinds of phenomena. First 
we have the true 3-body resonances with characteristic 3-body phase shifts [fTD] crossing 
7r/2, and with resonance behaviour in all partial waves. Secondly, we find structures 
(exhibiting fast growth of 3-body phase shifts up to 7r/2 in some partial waves) which 
are often caused by resonant and/or virtual states in binary subsystems. We call these 
structures 3-body virtual excitations. The analysis of the lowest partial components (those 
having the physically most transparent meaning) enables us in the HH method to obtain 
comprehensive insights into both the 3-body effects, and into the influence of resonances 
in binary subsystems on 3-body amplitudes. 

In [[19] we predicted and in |2(| explored in detail, in addition to the well-known 2 + 
resonance, (i) a second 2 + and a 1 + resonance in the 6 He continuum which both qualify 
as 3-body resonances, (ii) a soft dipole mode which does not but is a 3-body virtual 
excitation, and (iii) unnatural parity modes. Because of the halo structure of the ground 
state, there are peaks in the isoscalar responses of the soft monopole mode and soft 
dipole mode, even though there are no resonances in the low-energy continuum region. A 
higher-energy "breathing mode" appears as well, in the monopole continuum. 

Summarising the extended analysis of [pOfl , we show in Table | the positions and widths 
of possible resonances obtained by different methods. All of them give about the same 
positions, but different widths which should be testable experimentally. 

The way in which these structures could be experimentally observed depends on the re- 
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TABLE I. Comparison of resonance positions and widths of He. Results from the present 



Hyperspherical Harmonics Method (HH) and the Complex Scaling Method (CS) ]12,13| are 
shown, together with experimental data. Resonance positions are given relative to the theoretical 
ground state. 



action considered, as will now be demonstrated via nucleon charge-exchange and inelastic 
scattering to the 6 He continuum. A complicating feature is the overlapping of resonances 
and soft modes in the region of excitation energies 3-6 MeV, thus only a more detailed 
analysis of excitation functions and angular distributions may possibly distinguish those 
states. This will be a central issue of this paper, where we will show that even inclusive 
cross sections will be informative. 



III. MODEL DESCRIPTION 

In this section details are given on the physical ingredients of the model we have de- 
veloped for calculating inelastic and charge-exchange reactions to low-energy continuum 
states in 6 He. Structure and reaction scenarios are often intertwined in a very complicated 
way. In some situations, such as the very dilute matter of halo nuclei, the reaction dy- 
namics becomes simpler, and an approximate scheme using distorted waves is reasonable. 
In the DW framework, as follows from formulae below, the reaction amplitude has three 
ingredients: 

A. The structural information contained in the transition densities which describe the 
response of the nuclear system to an external field, 

B. The effective interactions between projectile and target nucleons, and 

C. The distorted waves describing the relative motion of projectile (ejectile) and target 
(residual) nucleus. 



A. Nuclear structure 



For description of the nuclear structure we have used the 3-body a + N + N model. 
In this model, the total wave function is represented as a product of wave functions 
describing the internal structure of the a-core and the relative motion of three interact- 
ing constituents (see appendix A). The method of hyperspherical harmonics (HH) (see 
Refs. pj, [IDLrB |) was used to solve the Schrodinger 3-body equation, for both bound and 
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continuum states. A modified SBB Gaussian type aN interaction |29|] with purely repul- 
sive s-wave component (Pauli core) fl[] and the "realistic" GPT NN interaction were 
used. 

We are now going to apply this model to continuum low-energy excitations above 
the (3-body) breakup threshold. The main model assumption about the factorization of 
the wave function into two parts suggests that low-lying nuclear transitions of interest are 
connected with excitations of the two valence particles in the halo outside the a-core. This 
assumption is physically reasonable for the low energy spectrum since a-core excitations 
must involve a significant energy transfer due to the particularly stable structure of the 
a-particle. 

The 6 He continuum reveals a variety of structures: the 1 + spin-flip resonance has an 
almost pure shell-model structure (P3/2P1/2)', the l\ and 2% resonances are on the other 
hand of strongly mixed nature; there are 3-body virtual excitation such as the soft dipole 
and monopole modes. 



For qualitative insight in the resonance structure, tables PHTTll give partial wave func 



tion norms in the interior region po < 15 fm for all resonances and the 1~ peak both in LS 
and (Jacobi) jj coupling scheme. These norms reflect the "eigen" resonance properties 
of any few-body system and measure the continuum strength accumulated in the strong 
interaction and centrifugal barrier regions. A hyperradial resonance wave function in the 
interior (p < p ) can be represented in a factorised form: 

x(r;E)~A ( E).*« (p ) = E _ (E C_. Tol2) .*« (p) (i) 

where \I/ R (p) has structure similar to that of a bound state. 

This general energy dependence is revealed by the reaction cross-section, but for wide 
resonances it is strongly influenced by the reaction mechanism. 
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TABLE II. Weights of the main components of interior parts of + wave functions of 6 He 



in LS and jj representation (ground state and O2 resonance at 5 MeV above the breakup 
threshold). 



The HH method is particularly suited for Borromean systems due to their simple 
asymptotic behavior. The physical characteristics of bound and low continuum states 
are concentrated in only a few wave function components corresponding to the lowest 
angular momenta and energy configurations of the 3-body system. This, combined with 
a convergence behaviour for ground state and resonances which is very much the same, 
preserves their relative position and enables us to avoid time consuming calculations. 
Thus we only take into account the hyperharmonics with hypermoments K < 6 that 
correspond to excitation energy up to about 10 MeV (up ~ K). Only the specific 3- 
body virtual nature of the soft dipole mode demands an substantially larger series of 
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TABLE III. Weights of the main components of the interior 2 + resonance state wave functions 
of 6 He in LS and jj representations at 0.8 and 3.0 MeV above the breakup threshold. 



hyperharmonics to achieve convergence. So we shall use the main components keeping 
in mind that the dipole mode will be somewhat shifted to lower energy within the peak 
width. 

B. Effective interactions between projectile and target nucleons 

The effective NN interaction V pt between projectile and target nucleons is a key point 
in the microscopic approaches to the description of one-step reactions. It differs from a 
free interaction since one of the nucleons is embedded in the nuclear medium, its motion 
being restricted by Pauli blocking and interactions with the nuclear environment. Usually 
these modifications are expressed by means of density dependence of the effective interac- 
tions, and a lot of work has gone into calculating effective interactions starting from the 
free one. As a rule, the calculations are based on nuclear matter with applications made 
to finite nuclei via local density approximations. These procedures include some uncer- 
tainties, which are especially troublesome when we deal with the lightest nuclei. Some 
physical situations exist, however, where the interaction dynamics simplifies, and simpler 
approaches can be used. At intermediate energies the impulse approximation has proven 
to be a very successful. In this situation the nucleon-nucleon collision energy is sufficiently 
large compared to binding energies, and the modification of the free interaction is not very 
significant. Using as effective interaction the free nucleon-nucleon t-matrix, that takes into 
account an infinite number of rescatterings between two nucleons interacting via a free 
NN potential, we obtain a complex, energy dependent interaction with parameters that 
can be extracted from analysis of experimental data on free NN scattering. 

Another simplified situation occurs in interactions with halo particles, because the halo 
particles have small binding energy and large probability to be outside the core of strongly 
bound nucleons. In the course of interaction with halo nucleons, small momentum and 
energy transfers are not blocked as is the case for interaction with core nucleons. As 
a result, the interaction with a halo nucleon is very similar to the interaction between 
two free nucleons, and we can use the free t-matrix interactions, in close analogy with 
the impulse approximation at intermediate energies. This approach is in the spirit of our 
model for nuclear structure, when the "active" part of the nuclear wave function is defined 
by the motion of halo particles. In concrete calculations of inelastic and charge-exchange, 
we use the i-matrix parametrization of Love & Franey with central, tensor and spin- 
orbit components. The contribution of an exchange knock-out amplitude is taken into 
account in the pseudopotential approximation p3| . 
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To be consistent with the 3-body a + N + N model of nuclear structure, we have 
also to take into account the aN-interaction between the projectile nucleon and a-core. 
For charge-exchange reactions at low excitation energy, only the halo nucleons are the 
active particles in our model. Charge-exchange with core nucleons must destroy the en- 
core and involve a large excitation energy, and was therefore neglected in our calculations. 
In inelastic scattering the aN-interaction can give a contribution also at low excitation 
energy. If the a-core is taken to have an infinite mass this contribution will be exactly 
zero, due to the orthogonality between initial bound and final continuum states of the 
3-body system. In real situations the center of mass of the total nucleus is somewhat 
shifted from the a-particle center of mass, and we have a finite contribution from this 
interaction. We expect that the role of direct aN-interaction increases with increasing 
transferred momentum. At this stage of our investigations, we neglect these contributions. 
Therefore, our calculations of inelastic scattering should be reasonable only for moderate 
transferred momentum. Physically this corresponds to the situation where the a-particle 
is a spectator, and experimentally it would be realized if only events with a-particles 
in forward direction were detected. It is necessary to underline that when we take into 
account the interaction with halo nucleons, the recoil effects are treated in an exact way 
because the wave functions which we use are defined in terms of the translational invariant 
Jacobi coordinates. 



C. Distorted waves 

For calculations of distorted waves we need to know the optical potentials describing 
the nucleon elastic scattering. For calculations of distorted waves we used a phenomeno- 
logical optical potential |35| describing proton elastic scattering from 6 Li at energy E p = 
49.5 MeV. The same optical potentials were used in both incident and exit channels. 



IV. REACTION FORMALISM 



Among the quasielastic reactions, the nucleon-nucleus inelastic scattering and charge- 
exchange reactions are the simplest and best understood. Experimental possibilities are 
now available for applying these reactions (in inverse kinematics on nucleon targets) to 
investigation of the structure of exotic halo nuclei. The cross section of quasielastic 
reactions, 

N + A -> N' + C + n x + n 2 

between a nucleon and a two-nucleon halo nucleus (core C, with A = C + n\ + n 2 in the 
g.s.) exciting the latter to the continuum, can be written in the form 



hvi 



J dk' N dk x dk % dk c S{E f - Ei)S{P f - PO | T fi | 2 (2) 



where Ei = e N + e A -, E f = e' N +e t +e 2 + e c + Q, Pi = k N + k A , P/ = k' N + ki + k 2 J t -k c 
are the total energies and momenta of all particles before and after collisions. In these 
expressions Q is the binding energy of the target nucleus in the case of inelastic scattering, 
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while it is the difference of binding energies of parent and daughter nuclei for a charge- 
exchange reaction. The relative incident velocity is Vi = — -, and //j = — — — is the 

Hi ' m N + M A 

reduced mass of the particles before collision. We will work in the center of mass (CM) 
coordinate frame (Pi = 0, k^ = -k A — hi), and use Jacobi coordinates for particles both 
in initial and final systems. The coordinates used are defined on Fig. [I] and are given by 





(a) 



(b) 



FIG. 1. Spatial coordinates (a) in nucleon-nucleus scattering, and particle momenta (b) in 
the two-neutron halo system. 
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In the CM frame Ef = e f + E K + Q, where £f = -, and \if 



m' N M A 



(3) 



is the 



reduced mass for the exit channel, while E K 



^x ~t~ ^1 



h 2 k 2 



+ 



m' N + M A 
is the excitation 



energy measured from the breakup threshold. Taking into account conservation of energy 
and momenta, the exclusive cross section (when energies and momenta of all particles 
are observed) is an average over initial and a sum over final spin orientation, and can be 
written as 



d 5 



a 



dQ x dQ y dQfde y dE K 



ft 4 h 



(2^^^2(^)IJe y (E K -e y ) 



1 



2(2J; + 1) 



fi 



(4) 



The factor y e y (E K — e y ) describes the distribution of energy between different modes 
of particle motion, and reflects the phase-space accessible for breakup. The matrix ele- 
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ment Tfi includes all the interaction dynamics, and is given in the Distorted Wave (DW) 
framework by 



T f i =< x { f ;li b (k f ), y { ~l m2 , mc (k x , ky) i e v pt i ^ JiMi , xtLiki) >, (5) 

t 

where Xif are distorted waves describing relative motion of colliding nuclei, ^j^a/, and 

^mi,m 2 ,m c ( , k y ) are the initial bound and final continuum nuclear states, respectively. 
Spin projections M a , Mi, Mb, m\, m<i and mc, together with relative momenta kij 
and k X) y characterize the asymptotic state of all particles taking part in the reaction. 
Continuum wave functions are a matrix in spin space, and contain the probability for 
spin-flip in the course of scattering. Since we are interested in studying the 6 He nucleus 
where the core is an a-particle with spin zero, for simplification of notation we omit 
everywhere below mentioning of the core spin projection mc- 

In eq. @ V pt is a local, effective nucleon-nucleon interaction between projectile and 
target nucleons p and t, expressed in terms of central, spin-orbit and tensor components, 



V(r pt ,p pt ) = £ |E^tM^ • of + t T LS ( r Pt)L ■ S + tl{r pt )S pt {r pt ) U T • 7? 

TVS ) 



(6) 



where r pt = r p — r t and p pt = \{p P — Pt) are relative distance and momentum (wave 
number) between two nucleons, S = \(<r p + cr t ) and L = r pt x p pt are operators of total 
spin and orbital angular momentum of the two nucleons, p L = — zVi, erf = 1, <Tj for S=0, 
1 respectively. S pt (r pt ) is the tensor operator 

q f4t s _ H(Tp-r pt )((Tfr pt ) 

r pt 

Using the method of hyperspherical harmonics, the nuclear wave function above the 
breakup threshold can be written as follows (for details see Appendix A): 

*™ = E < s i™i « 2 m 2 I S f M Sf >< L f M Lf S f M Sf | J f M f > (8) 

■y,J f ,M f ,M Lf 

where 7 is an abbreviation for a set of quantum numbers 7 = {Kf,Lf,Sf,l x ,l y }, which 
characterizes the relative motion of the three particles flying apart. The continuum wave 
function ^f^j-Mf V, depends on the quantum numbers 7, Jacobian space coordinates 
(x, y), nuclear excitation energy E K (expressed by the hypermomentum k), and the total 
angular momentum Jf and its projections Mf, 



%,j f , Mf (x, y, k) = j—^ JlXKn>j v f ,K f i x i y M Tjfhvs'Mf&s) ( 9 ) 
The transition amplitude Tfi can be further decomposed according to eq. (||), 
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T fl = < Sim lS2 m 2 | S f M 8f > < L f M Lf S f M Sf \ J f M f > * Mhf T,(k x , k y 



1 ,M Lf ,M Sf 



(10) 



where T^kx, k y , k) formally now has the same structure as any two-body amplitude for 
excitation of a nuclear state with total momentum Jf,Mf, excitation energy E K , and a 
fixed state of relative motion of breakup fragments defined by the quantum numbers 7: 



T^k x ,k y ,K) = < XfMM'^-rJf'MtiK) I Y. V Pt I ^JiMi,Xi!M a ( k i) > 

t 

E < JM jm I J f M f > T^ Mb ' m {k x , k y , k) (11) 



To calculate the reaction amplitude we use the partial wave decomposition for the dis- 
torted waves XiMai^i describing the relative motion of the projectile nucleon and the 
center of mass of the target nucleus: 

= 7^~E < Lm la S a M a \ 3a m a > i l <* Y* mi (k { ) Xi aja (h,r p ) \j a m a ) (12) 

where \j a m a ) = E < Lm' la S a M' a \j a m a > Y lam ^(r p ) \S a M' a ), and \S a M' a ) is the projectile 
spin function. Nuclear formfactors can be defined as follows 

< j b m b , J f M f I J2 V Pt I JiMi,j a m a > = E < J i M d m I J f M f > ( 13 ) 

t Isj 



I _]\l+s-j+rn a 
-I S Zl / I W.(& \ ft* . L' m \ P is J C K r _ 



8 < I [Yi(r p ) <g> a p ] jm I j a m a > F/j a (K, r p , — ) (14) 



where s = V2s + 1. In the case of nuclear excitations of normal and unnatural parity 
the explicit formulas for radial formf actor are given in plj and Appendix B. Taking 



JbJa 

nM a ,M b ,rn , 



into account these definitions the reaction amplitude T f" b ' m (k x , k y , k) can be written 
in usual form 



1 1,3 {Kx,Ky,K)- 2^ 1 1 l a m la *l b m lb K K f ) 1 l a ja,hjb IJ^ 



h S b j b 

x < l a mi a S a M a \j a m a > < l b m h S h M h \j b m h > < j b m b jm\j a m a > < l b 0l0\l a > {I s j 

la S a j a 



where the radial integrals Ii*j a ,i b j b are defined as follows 

1 d I 

- Xi bjb (kf,r p ) F/ a J jb (K,r p ,—) - 



is- f°° 1 ■ (9 1 

^ZaMb = J Q dr P r l - Xi bjb (kf, r p ) F/J jb (n, r p , — ) - Xi aja (k h r p ) (16) 



It is useful to compare the expression (|i~0|) for the breakup amplitude T/j with the 
amplitude for a usual two-body reaction which has the same structure as amplitude flllD - 
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In the 3-body case, the amplitude Tfi has additional degrees of freedom which are man- 
ifested as dependence on angles Q§ = {a,k x ,k y }, where sin 2 a = — In contrast to a 

E K 

traditional two-body approach, the exclusive cross section (proportional to | Tfi | 2 ) con- 
tains an incoherent sum over total spin Sf but a coherent sum over total transferred j 
and final J/. Consequently, we expect that the exclusive cross section will be especially 
sensitive to the correlations in the nuclear structure. 

The different exclusive cross sections and correlation distributions will be considered 
elsewhere, here we restrict ourselves to inclusive cross sections. To calculate the double- 
differential inclusive cross sections when experiments measure the energy and angle for 
one particle, we must integrate the fivefold exclusive cross section (f|) over the unobserved 
coordinates of breakup particles (angles x = Q x , y = Q y ), and over various distributions 
of relative energy e y between fragments: 

— — — - = / de v / dVt x dVt v -— — — — — — - — — — (17) 

dilfdE K Jo J dil x dYL y d\Lfde y dE K 

Using the following orthogonality properties of the hyperspherical harmonics, 

rE "de y y/e v (E K -e v ) J dQ^yf^ (^jft^) 

= ^ElS K r K Sit L 5 M ' LML 5i' x i x S Vy i y , (18) 

the inclusive cross section now becomes an incoherent sum over total transferred j and 
final Jf angular momenta, and different 7-components of the final target state are excited 
independently of each other. Thus 

d?a A HiHfkf x ^ (2J/ + 1) 

" ( 27r ) -T^IT £ 



dQ f dE K K ' h A h ij ^ 7 (2J i + l)(2j + l 

1 \ - I rpm a m b m (1 , x , 2 . j-ai^xi^y 



E ir- m Hfc/,fc^)i 2 4^(^)i. (19) 

In this expression the factor E 2 , which originates from the 3-body phase volume, guaran- 
tees the correct cross section behavior at the breakup threshold. From ( |T9"D it also follows 
that, due to the averaging procedure, we lose information about correlations in relative 
motion of the breakup particles (which were defined by ^^(^5) hyperharmonics); 
remnants of the complex dynamics that governs the particles motion are kept only in dif- 
ferent shapes and strengths with which various components of final states are distributed 
over excitation energies. One may hope that in the differential inclusive cross sections, 
due to specifics of reaction mechanisms, we can under certain conditions enhance the 
excitations of some of the components and thus still obtain valuable information about 
structures of halo nuclei. 

The inclusive excitation cross section can be obtained by integrating over all ejectile 
angles Qf~. 

da f d 2 a 

dE K = J dnf dnJdE K (20) 

This cross section describes the distribution of total strength of different excitation modes 
over energy spectra in quasielastic reactions. 
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V. RESULTS 



With the model described above, we have calculated the excitation (eq. ^Dp and 
double-differential (eq. |19[) inclusive spectra in the CM system for the charge-exchange 
reaction 6 Li(n,p) 6 He and the inelastic scattering 6 He(p,p') 6 He at E N = 50 MeV, with 
excitation of different Jj 1 = ± ,1 ± ,2 ± ,3 ± low-energy states of 6 He. The corresponding 
cross sections are shown in Figs. In the figures E* is the nuclear excitation energy 
measured from 6 He ground state. For inelastic scattering the initial target state of 6 He 
has Ji = and the total j transferred has a unique value and coincides with Jj of the 
final state. For charge-exchange on 6 Li the situation is more complex: since J, = 1 for the 
6 Li ground state it is possible to excite final states of 6 He with definite J/ by different j 
transfers. All values of j allowed by angular momentum conservation (Jj — J j = j) were 
taken into account in our calculations. It follows from equation ( |19"D that every j gives an 
independent contribution to the inclusive cross sections. 

Our main goal is to demonstrate, that even in the simplest inclusive experiments 
it is still possible to extract information about structures in the continuum by detailed 
examination of both excitation and differential cross-sections. 



A. Two test cases for the model 

Two cases were used to check the model and consistency of our reaction continuum 
calculations. The sharp 2f resonance at 1.8 MeV was used in the first. This resonance 
resembles a usual bound state and can be described with good accuracy by calculating it 
with a boundary condition under the barrier corresponding to a discrete state. We use 

this to calculate the differential cross section — . Next we calculate the double-differential 

d£l 

d 2 a 

— — - — cross sections for the 2f resonance at different E* and after that integrate over 
dD.dE* 1 

E* across the resonance. In fact, calculating a resonance width and cross section at 
peak position, energy integration has been done analytically since the resonance has the 

Breit-Wigner form (we checked it). In both calculations we got the same results for — -. 

The second way is to compare our calculation with known experimental data for 
excitation to the continuum. In work [34| the reaction 6 Li(n,p) 6 He at neutron energy 118 



MeV was measured. The proton energy resolution in the experiment was ~ 2.3 MeV. 
This should be kept in mind when comparing with the reported differential cross sections 
for transitions from the 1 + ground state of 6 Li to the + ground state and 2f resonance 
(1.8 MeV) of 6 He. Fig. shows the corresponding experimental data plotted together 



with our calculations using the Love and Franey t-matrix interaction |32] at 100 MeV and 
an optical potential describing proton elastic scattering from 6 Li at 144 MeV. A good 
description for the shape and absolute value of the differential cross section to ground state 
was obtained and also a reasonable agreement with the data on the 2 + resonance. The 
2 + angular distribution has a characteristic form corresponding to a transition of mixed 
angular momenta. To demonstrate this we show in Fig. |2|b the separate contributions 
from transitions with total transferred j equal 1, 2 and 3 by dashed, dotted and dashed- 
dotted lines, respectively. The transition with j = 1 includes transfer of relative orbital 
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momentum I = and determines the cross section at small angles, the others with j = 2 
and 3 have / = 2 and dominate at larger angles. 

Thus the reliability of our approach was confirmed, arguing for our predictions for 
low-energy excitation spectra, for which the model was developed. 




FIG. 2. Angular distributions for 6 Li(n,p) 6 He at 118 MeV. The dashed, dotted and 
dashed-dotted lines show the contributions from j = 1,2 and 3, respectively. The experimental 
data are from ref. [34|. 



B. Inclusive excitation spectra 

1. Partial content 

The inclusive excitation spectra (Fig. |3], thick solid line) for charge-exchange and 
inelastic scattering reveal two distinguished bumps in the low-energy total spectrum : The 
first, narrow at excitation energy ~ 1.8 MeV and the second, broad at ~ 4.5 MeV. To 
understand the nature of these structures, the left side of Fig. ^| shows the decomposition 
of the total spectra into contributions from excitations of different partial components 
JJ of 6 He. The 1~, 2 + and + excitations are given by thin solid, dashed and dotted 
lines, respectively. The dot-dashed line shows the contribution from 1 + excitation for 
charge-exchange and 3~ for inelastic scattering. Contributions from other partial waves 
are less significant and not given in the figure. 
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FIG. 3. Multipole Jj f and spin decomposition (left and right sides) of inclusive proton energy 
spectra from 6 Li(n,p) 6 He* (top row) and 6 He(p,p') 6 He* (bottom row) reactions. 

The first narrow peak is the well-known 2f resonance in 6 He. The broad bump has 
a more complex structure. A mixture of different excitations is responsible for the total 
shape; a second 2 + resonance and concentration of low lying strength of l - and + 
excitations dominate the spectrum. The double-hump shape of 2 + excitations is the 
most remarkable feature of the low-energy spectrum. The strength concentration of 1~ 
transitions at E x ~ 4 MeV is the other peculiarity. The behavior of other excitations, for 
example 3~, is different. It smoothly increases from threshold, and in the case of inelastic 
scattering gives a significant contribution at higher excitation energy. 

The excitation spectra for both reactions have qualitatively the same gross structure, 
but the absolute cross sections are a few times larger for inelastic scattering than for 
charge-exchange. 

Nuclear reactions in which halo nuclei take part serve, due to somewhat different dy- 
namics, like filters and could lead to different multipole composition in observed excitation 
structures. In inelastic scattering the dipole mode dominates while in charge-exchange 
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the 2 + resonance is about 50% larger. 

The pronounced 1~ nuclear excitation has similarities with electromagnetic response 
for the soft dipole mode, prevalent in Coulomb breakup on heavy targets. Fig. [| shows 
theoretical cross sections for Coulomb breakup (dotted line) of 6 He on gold at 63 MeV/A 
and inelastic proton scattering (solid line) with 1~ excitations, arbitrarily normalized. A 
semiclassical description is used for the Coulomb dissociation process and our model for 
the electromagnetic dipole response in 6 He 0. Both processes show strength accumula- 
tion in the same energy region and hence, we should expect no matter which excitation 
mechanism dominates, a similar behaviour for excitation functions. In an elegant ex- 
periment on 6 He breakup reaction at 63.2 MeV/nucleon on Al and Au targets [24] with 
registration of 7— rays in coincidence, similar behaviour of a— particle distributions was 
found for both targets. For the light Al target, nuclear mechanism is believed to give the 
main contribution to the spectra while for Au the EM dominates. Our theoretical results 
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explain qualitatively the observed similarity. 



2. Spin structure 

The composition of spectra, or relative role of excitations of various Jj f , is as discussed 
above different for the two reactions. In charge-exchange a relatively larger number of 
states was excited with about equal intensity, while inelastic scattering is more selective. 
To better illustrate this point, the right side of Fig. |3| shows separately the contributions 
from excitations of 6 He states with total spin Sf = (dashed line) and 1 (dotted line). For 
inelastic scattering the excitations with Sf = dominate the spectrum, while for charge- 
exchange both contributions become comparable. This is a reflection of specific reaction 
mechanisms. In inelastic scattering the S = 0, T = component of effective interactions 
is the biggest one, while in charge-exchange only isovector components play a role and in 
the charge channel the effective forces with S = and S = 1 are comparable in strength. 
The relative role of different components of effective forces depends on collision energy 
and so the ratio between excitations of the various structures will change accordingly. 

C. Double-differential inclusive cross sections 

1. Fixed angle 

Excitation functions, measured at a fixed angle, can serve as a filter for selecting partial 
waves with definite multipolarity and therefore make it possible to extract information 
on resonances in complex situations like that described above. 

Fig. |5| shows spectra for charge-exchange at different exit proton angles. The total, 
1~, 2 + , + and 1 + spectra are denoted by thick solid, thin solid, dashed, dotted and 
dashed-dotted lines, respectively. The double-hump shape appears at all angles, but 
the composition of the second bump depends on scattering angle or angular momenta 
transferred in the reaction. At 0° the excitation of + , 1 + and 2 + dominates the spectrum. 
All are populated by strong transitions with relative orbital momentum I = 0. With 
increasing angle the excitation of 1~ grows and at 10° all these states are important. 
At 20° the 1~ and 2 + are most pronounced. At 40° the absolute cross section has been 
reduced to half value, and all energy spectra except 2 + have flat distributions. It is also 
interesting to note that for 1~ excitation the main contribution comes from the transition 
with j — 2, s — 1 and 1 = 1. The dominance of spin-flip transition is due to the structure 
of initial and final states where components with Si = 1 and Sf = prevail. 

Fig. ^| shows the analogous spectra for inelastic scattering. All lines mean the same 
as in the previous figure, except that the dotted line denotes 3~ excitation for 9 = 20° 
and 40°. At 0° the + is excited very effectively in the region of the second bump. With 
increasing scattering angle 1~ becomes pronounced. It dominates the total spectra at 10° 
and 20°. Contributions from 2 + and 1 + additionally increase the width of this bump. 
At 40° the total cross section is again diminished, with 2 + and 1~ excitations being the 
largest. At higher excitation energy 3~ now gives a significant contribution to the total 
spectrum. 
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From the results represented in the figures, it follows that the second bump structure 
in the low-energy part of the spectra is a complex mixture of various excitations of the 
6 He nucleus. The way these excitations are revealed depends on the external fields (or 
reactions) applied to the system. This is a characteristic feature of continuum excitations 
without sharp resonances. 



2. Fixed excitation energy 

We discussed above the dependence of the differential cross sections on excitation 
energy for fixed scattering angle or momentum transfer. It is also interesting to compare 
the behavior of the cross sections at fixed excitation energy for different scattering angles. 
Fig. [7] shows angular distributions for the (n,p) reaction, for a few excitation energies that 
cover both sides of the second bump. The values of the transferred momenta q=ki — kf 
(in units of fm _1 ) are shown on the top abscissa, corresponding to the scattering angles 
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6 He(p,p') 6 He* 




E* (MeV) E* (MeV) 

FIG. 6. Inclusive differential proton energy spectra from 6 He(p,p') 6 He* for four values of the 
scattering angle 6 cm . See the text for further details. 

shown at the bottom abscissa. The thick solid line shows the total cross section. On 
the left side of the bump (E* = 3.5 MeV, Fig. 0a) the differential cross section has an 
asymmetric bell shape with maximum at about 15°. Going to higher excitation energy 
(E* = 4.1 MeV, Fig. ^b) through the bump maximum(E* = 4.5 MeV, Fig. [7|c) to the 
right side (E* = 5.2 MeV, Fig. |7|d), the cross section shape changes smoothly and becomes 
gradually more flat with a plateau from 0° to 20° on the high energy side. This shape 
modification becomes transparent if we examine how the contributions from excitations of 
different Jj f change with energy. The thin solid, dotted, dashed and dash-dotted curves 
show the contributions of + , 1~ , 1 + and 2 + , respectively. We see that at small angles, the 
dominating states are all those (0 + , 1 + and 2 + ) which can be reached by transitions with 
zero relative orbital momentum. For a range of somewhat larger angles the contribution 
of 1~ is most significant. The + and 1 + have smoothly falling angular distributions, 2 + is 
more flat due to the already mentioned mixing of transitions with different j. Hence, the 
interplay between + , 1 + and 2 + excitations, which together create a smooth background 
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FIG. 7. Proton energy spectra from the 6 Li(n,p) 6 He* versus scattering angle 8 cm for four 
excitation energies E*. 



with highest cross section at small angles, and the 1~ peaking at 20°, define the total 
shape. The competition between them is responsible for the modification of this shape 
with excitation energy. As a result, we get a flat total distribution extending over a rather 
wide angular range on the high-energy slope of the second bump. These results are in 
qualitative agreement with experimental data on the 6 Li( 7 Li, 7 I?e) 6 He* reaction [27,2S| if 
we scale angular distribution according to the transferred momentum q. 

The corresponding data for inelastic scattering are shown in Fig. || The thick solid, 
thin solid, dashed and dotted lines again show total, + , 1~ and 2 + cross sections, re- 
spectively. For inelastic scattering, in contrast to charge-exchange, the total cross section 
remains bell-shaped at all excitation energies of the second bump. This is caused by the 
dipole excitation which dominates the spectra. The contribution from + is also signifi- 
cant, especially at small scattering angles, counteracting the drop of the total cross section. 
The excitation of 2 + in (p,p') does not play the prominent role it does in charge-exchange. 
Together, + and 2 + , create a smooth background in angular distributions. These dif- 



19 



6 He(p,p') 6 He* 



0,11 0,25 0,45 0,67 0,11 0,25 0,45 0,67 0,87 



> 

B 
E 



C! 

■a 

*UJ 

«? 

T3 



1,5 



1,0 



0,5 



0,0 



°» 1,5 
> 
a> 
E 



£3 
§ 

■a 

*UJ 

«? 



1,0 



0,5 



0,0 



1 1 1 1 1 1 1 

E* = 3.5 MeV 


— « — i — «■ 
E* 


1 ■ r 

= 4.1 MeV 




(a) 






(b) 


/ r \ 




r \ 

X \ 

\ 

\ 

\ 

\ 

\ 

\ 




/ / N \ 
/ ' V \ 

/I \ \ tot 

i \ \ 

/ \ \ ■ 


~^^° + 


tot 

\ 


/ ^ N ^ 


T 






-|- | 1 1 1 1 1 

E* = 4.5 MeV 

(c) 


1 1 >■ 

E* 


— i — i — y 
= 5.2 MeV 


(d) 


/ \ tot " 










/ 

/ 

/ 


r 

s. 

+ 


\tot 

N. 









' 10 20 30 10 20 30 40 



FIG. 8. Proton energy spectra from 6 He(p,p') 6 He* versus scattering angle 6 cm for four exci- 
tation energies E* . 



ferences in the two reactions are due to two reasons: i) because of different structure of 
initial states we need different operators to excite the same final state in 6 He and ii) in 
addition to isovector forces in charge-exchange the strong isoscalar NN interaction acts 
between target and projectile nucleons in inelastic scattering. 

To further understand the nature of continuum excitations it is also useful to make a 
comparison between angular distributions for the two 2 + resonances: the first one being 
narrow and the second broad. Fig. |^a shows the differential cross sections for 2f at three 
excitation energies: approximately at peak (solid line) and at energies shifted from the 
peak position a half width to the left (dashed line) and to the right (dotted line). We 
see that all angular distributions have identical shape through the resonance. Fig. |9|b 
shows separately the contributions from excitation of the three main components of the 
6 He wave function ( 1 (dashed) is for L = 2, S = 0, l x = 0, l y = 2; 2 (dotted) is for 
L = 2, S = 0, l x = 2, l y = and 3 (dot-dashed) is for L — 1, S — 1, l x — 1, l y — 1) to 
the total (solid line) 2^ cross section at peak position. Again, the angular distributions 
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FIG. 9. Proton energy spectra from 6 He(p,p') 6 He* for the first narrow 2 + resonance: (a) The 
cross sections at peak position and at energies shifted from the peak by a half-width to either 
side, (b) Partial contents of the peak cross section. For details, see the text. 



for all components have the same shape. For the broad resonance the picture is different, 
as shown in Fig. 10a where solid, dashed and dot-dashed lines show the total angular 
distributions for l\ at peak position and shifted from it by a half-width to the left and 
right, respectively. For the broad resonance, the shape of the differential distribution is 
changed through the resonance. Fig. |lO"lb,c,d shows decomposition of the total distribution 
into contributions from main components. The notation is the same as in Fig. ^jp. For 
2% the main contribution comes from excitation of final quantum numbers L = 1 and 
S = 1 (curve 3). The shape of this component is only slightly changed when going 
across resonance: the maximum shifts by about 2° and the width becomes narrower on 
the high-energy side. The shapes of other components experience dramatic changes: the 
interference pattern in the angular distributions has a different character on opposite 
sides of the resonance. Usually the resonance amplitude, as function of energy, can be 
separated into a smooth background and a resonance part. It is reasonable to assume 
that for a sharp resonance the background part remains more or less constant over the 
resonance width and all energy dependence is only in the resonance part. As a result, the 
shape of angular distributions does not change over the resonance. For broad resonances 
the background part may change, and interference with the resonance part can produce 
different angular distributions. For dominant components the role of the background part 
is relatively small, hence shape variations are not pronounced. For smaller components 
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FIG. 10. Proton energy spectra from 6 He(p,p') 6 He* for the second 2 + resonance: (a) The 
cross sections at peak position and at energies shifted from the peak by a half-width to either 
side. The partial contents of the corresponding cross sections are shown in (b), (c) and (d), 
respectively. See the text for further details. 

both parts of the amplitude are comparable, and can give different angular distributions 
on opposite sides of the resonance. 



D. Transition densities to bound and continuum states 

Another interesting illustration can be obtained from comparison of transition densi- 



ties to bound and continuum states. As an example, Fig. |ll] shows transition densities in 
momentum space from the 1 + ground state of 6 Li to ground and continuum + states of 
6 He for transferred orbital, spin and total angular momenta equal 0, 1 and 1, respectively. 
The continuum energy was chosen as 6 MeV, where excitation of the continuum + in 
charge-exchange is largest. Since transition densities to continuum are complex, we show 
only absolute values. Curve 1 is for transition to ground state, curves 2,3 and 4 are com- 
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ponents of continuum 0^ with quantum numbers (K — 2, L — 1, S — 1, l x — 1, l y — 1), 
(K = 2, L = 0, S = 0, l x = 0, l y = 0) and (K = 0, L = 0, 5 = 0, l x = 0, l y = 0), respec- 
tively. We see that between bound states the transition density has a unique spectral com- 




FIG. 11. Absolute values of the charge-exchange momentum space transition densities 
pl—Q S= i j—i(q) to the + states of 6 He are shown. Curve 1 is the 0^" density, curves 2-4 are 
the O2" density components. For details, see the text. 

position, which coincides with that known from the electron scattering Ml form-factor. 
For transitions to the continuum, there are various components with different spectral 
forms. In the bound transition, similar components with the same quantum numbers 
also exist, but since the bound state presents a unique structure all these components are 
organized in a unique way, and give a joint system response to external perturbation. In 
the continuum the various components correspond to different modes of relative motion 
between breakup fragments and should, in principle, be accessible to measurement. They 
will be excited differently by different reactions and the response will depend from the 
external fields applied to the system. Only in the case of sharp resonances (which in many 
respects are similar to bound states and represent to a large extent the internal property 
of the system) will the response be more or less the same. 



VI. CONCLUSION 

We have developed a four-body distorted wave theory which is appropriate for analysis 
of nucleon-nucleus reactions leading to continuum excitations of two-neutron Borromean 
halo nuclei. Spatial granularity of the halo bound state and the final state interaction 
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in the 3-body continuum was fully taken into account by the method of hyperspherical 
harmonics. The weak binding and dilute matter of halo systems enabled us to use a 
free NN t-matrix for the interaction with halo nucleons. Although applicable to any 
two-neutron Borromean halo nucleus, the A = 6 nuclei were again chosen as benchmark 
systems. For these nuclei we have the most complete knowledge of the binary subsystems. 
Experimental investigations of these nuclei are also currently performed or planned. As 
an initial check the model was successfully tested against data for elastic 6 Li(p,p) 6 Li, 
inelastic 6 Li(p,p') 6 Li (0 + , 3.56 MeV) and 6 Li(n,p) 6 He (gs. and 2f, 1.8 MeV), which has 
been available for a decade. 

A detailed study of inclusive excitation and differential cross sections for inelastic 
6 He(p,p') 6 He* and charge-exchange 6 Li(n,p) 6 He* reactions at beam energy 50 MeV was 
performed. The theoretical low-energy spectra exhibit two resonance-like structures. The 
first (narrow) is the excitation of the well-known 2± resonance. The second (broad) bump 
is a structural composition of overlapping soft modes of multipolarities l~,2 + ,l + ,0 + 
whose relative weights depend on transferred momentum and reaction type. Recent ex- 
perimental data on heavy- ion charge-exchange reactions [p7[p8f , although sparse, confirm 
the existence of the second structure. 

The soft excitations of different multipolarities have a concentration in a relatively 
narrow energy region near the 2f resonance. This poses a challenge. Nuclear reactions in 
which halo nuclei take part serve however, due to differences in reaction mechanisms, like 
filters emphasizing different multipole components in the observed excitation structures. 
To some extent we may exploit this to our advantage. 

Thus comparison of (n,p) and (p,p') shows that the excitation cross section for inelastic 
scattering preferentially selects the 1~ component. Hence (p,p') is the most promising tool 
for studying the soft dipole excitation mode. 

Double differential distributions for the broad structure show that association of the 
observed structure with excitation of a unique multipolarity would be misleading. This 
is especially so for charge-exchange 6 Li(n,p) 6 He*, where a flat shape of the total angu- 
lar distribution extending outside forward angles, is due to mixing of excitations with 
different multipolarities. Under favorable conditions, measurement of spectra at definite 
momentum transfer makes it possible to extract information on individual resonances in 
complex situations like the one described above. 

Our results on charge-exchange are in qualitative agreement with experimental data 
on the 6 Li( 7 Li, 7 Be) 6 He* reaction [271 28 1 if we scale angular distributions according to the 
transferred momentum. Forward angles are most important for partial analysis, but in 
both experiments there is not enough statistics in this region for more definite conclusions 
on the resonant structure of 6 He continuum. Since all resonant states are concentrated 
in the vicinity of the extremely pronounced 2f state, high resolution experiments with 
detailed angular distributions will be needed. 

The model we have developed allows us to calculate cross sections for kinematically 
complete experiments when characteristics of four particles are measured. Hence we can 
study a variety of correlations existing in Borromean halo nuclei that could not be seen 
in the inclusive observables. An analysis of different exclusive cross sections of nucleon- 
nucleus reactions with excitation of the 3-body continuum of 6 He is in progress. 
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APPENDIX A: 



Within the cluster representation (for details see Refs. m FOlJlB"! ), 3-body bound and 
continuum state wave functions (WF) have the product form 

| $ > = exp(zK • R) $ c (Cc)*jm ( A1 ) 

where $>c{Cc) is an intrinsic core WF, while ^ T JM is the "active" part of the 3-body WF 
carrying the total angular moment J, its projection M and total isospin T. This part 
depends on relative coordinates and cluster spins (suppressed in our notations) and it is 
the object of the calculation. K and R are momentum and coordinate of the center of 
mass of the nucleus A, respectively, 

Translationally invariant normalized sets of Jacobi coordinates X3 and y 3 are defined 
as follows 



x 3 = \/A 12 (r 2 - r X/ 

Air 1 + A 2 r 2 



V3 = sjA m c (r c - Ai + A2 ), (A2) 
R = - (Ain + A 2 r 2 + A c r c ). 

Here A X2 = AiA 2 /(Ai + A 2 ) is the reduced mass of the (12) subsystem in units of the 
nucleon mass m^, ^4(i2)c = (Ai + A 2 )Ac/ (Ai + A 2 + Ac) is the reduced mass of the (12) 
cluster with respect to the core C, and A = Ai + A 2 + Ac- Notice, that is co-linear 
with rc — R- Alternative sets (x-i,yi) and (#2, 2/2) of Jacobi coordinates are obtained 
by cyclic permutations of (1,2,C). The set of Jacobi momenta ^3, P3 and K conjugate to 
x 3 , y 3 and R is defined by the relations, 

k c ki + k 2 



K = ki + k 2 + k c , 

where ki, i — 1, 2, C are the particle wave numbers in an arbitrary frame. The Jacobi 
momenta q%, p 3 are connected to k x and k y defined in (|3|) with simple relation 
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k x = \[Ai2 <?3, 

k y = JA (12)C Ps, (A4) 



We use hyperspherical coordinates p, a, 9 X , <p x , 9 y , <p y , where (9 X , (fi x ) and (9 y , <p y ) are 
angles associated with the unit vectors x and y, and 

p = (x 2 + y 2 ) 1 ' 2 = ( ]T Mr, - R) 2 ) 1 / 2 , a = arctan(x/y). (A5) 

i=l,2,C 

The collective variables a and p are called hyperangle and hyperradius. The last variable 
is rotationally and permutationally invariant, having the character of total moment of 
inertia or a weighted measure of distances in the 3-body system. The corresponding 
conjugated momenta are 



K 



= (q 2 +P 2 Y /2 = h-\2m N I E K |) 1/2 , a K = arctan(g/p) (A6) 



where E K is the total 3-body energy. Since we introduce a new degree of freedom its 
corresponding conjugated quantum operator has eigenvalues K = 2n + l x + l y called 

the hypermoments. Hyperspherical harmonics (HH) lp^ v {a) ■ Yi xmx (Q x ) ■ Yi yrny {Vt y ) are 
eigenfunctions of of this operator. 

We seek our bound-state and continuum wave functions in the form of expansions on 
a generalized angle-spin basis (LS coupling) 

tSw^) = \y l K l m>)®xs\ JMj (A7) 



with HH defined as 



= fP(a) Y lx (Q x ) ® Yi y (Q y ) lm (A8) 



Here the a, 9 X , <p x , 9 y and (fi y variables are denoted collectively by Jl 5 , X s is a spin 
function, ...(g)... means vector coupling, 

Ai <g> B x ] = J2 < lm i | jm > A lmi B Xmx (A9) 

The relative orbital momenta l x , l y , couple to the total orbital momentum L and its 
projection M. Hyperangular part of HH has the following explicit form 

l/) l ^(a) = N 1 ^ (sin a)'- (cosa)'" P£gffi +1/2 (cos 2a), (A10) 

2 

where -P" ,/3 are Jacobi polynomials and N l ^ ly is a normalization factor. 
For bound states the internal WF in LS coupling has the form 

*Im = ^ E X£L>) T^ i5Mj (^) A TMt . (All) 
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where 7 is an abbreviation for a set of quantum numbers 7 = {X, L, S, l x , l y }. For 
continuum states we have the following form 
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/ 



^ ' JM ~ ( K p) 5 / 2 ^ XKl' x l y ,K'l' x l' y ( K P) ^JKLSMji^) 

x < L'M' L S'M' S | JM > Xtm t (A12) 

with normalization condition 

^J^dxdy = «T 5 <5(/t' - - fig) = <S(g' - q)6(p' - p) (A13) 

The WF \I/J M is a solution of the 3-body Schrodinger equation 

(f + V- E)^ T JM = , V = V 12 + V 1C + V 2C , (A14) 

where V^- is the interaction potential between particles i and j. After separating out the 
hyperangular parts of the WF we obtain a set of coupled equations similar to those for a 
particle moving in a deformed mean field. 

In case of neutral particles the bound hyperradial WF for Borromean nuclei has a true 
3-body asymptotics 

X> - 0) ~ p K+5/2 ; X> - oo) ~ exp(- K p) (A15) 

For continuum WF the boundary condition at the origin coincides with that for the bound 
state, while for chargeless particles at p — > oo it is 

X 7) yM ~ P 1,2 {H K+2 {K P )5 ia , - S^H+ +2 (k P )) (A16) 

Here if "and are Hankel functions of integer index (n — K + 2) with asymptotic 

~ — — exp(=piKp), describing the in- and outgoing 3-body spherical waves, S^y is the 
v P 

^-matrix for the 3 — > 3 scattering. 

Wave functions discussed above are characterized by the total angular momentum J 
and its projection M. Due to rotational invariance the continuum wave functions with 
different J are dynamically decoupled and can be calculated separately. For transition 
densities we need the 3-particle scattering states m2 (k x , k y ) in other representation 
characterized by k x and k y momenta of relative motions and projections m\ and m 2 of 
particle spins on a chosen direction. They can be written as follows 

^Li,m 2 = / N5/2 E XKl x l y ,K'l' x l' ( K P) f JKLSm{^) 3 / 'k'L> M'JS^T) 

\Kf>) ,J,M,M' L 

x < L'M' L S'M' S \JM>< s 2 m 2 | S'M' S > X TMt . (A17) 

The transition density describes the system response to a zero-range perturbation, and 
can be expressed as a matrix element between initial bound and final continuum states 

f%£ T = < I E ^^[W ® at] rZ T (t) | J,M 4 > 

t=l,2 r * J 

= -i- < JMi jm | J/M/ > E < s i m i s 2"i2 | S' f M' Sf > 

j f> M f J f y P M ' Lf ' M s f 

x < L'jM' Lf S>M' Sf | JfMf > y°(Z/*, f m p l *F(r, «) (A18) 
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The easiest way to calculate space integrals in p l y' T (r, k) is to do it in a coordinate system 
where radius r t is colinear with the y-coordinate: 



/ Ai + Aj 



r t = a k Vk = \j A ^ Vk (A19) 

In our case we must rotate from the initial Jacobi coordinate system (basis x 3 ,y 3 
with A 1 = A 2 = 1) to the alternative similar sets {x^y-i) or (x 2 ,y2)- Hyperharmonics 
y l KLM transform under this rotation through Raynal-Revai coefficients 



y l KL V M{^) = £ < C l'y\lx, K >KL ( A2 0) 



/' /' 



Using the following definition of reduced matrix elements 

i A i -r < ji m i l m \ 1f m f > ■ i i A ii- / , 

<3fm f \O jm \J i m i >= J Ji } < J/HOjHfl > (A21) 

Jf 

and with the necessary summation over Clebsh-Gordon coefficients, the radial part of 
transition density matrix elements p l y' T (r, k) is 

p l y ,T (r,K) = £ < l x f l y f \ l Xfly f >K f L f < I'x/y^xjyi > K t L, 

-vj- -v- 7" I" I" I" 

If ! Ht l x f >'j/y I'j/j 

x < ^11^(1)11^ >< > <% r ^ 

f 1" 1" T ~\ I ^ ^ 

x *W l xi p J L . L/ l I (A22) 



Si Sf s 



> 



The factor (— l) s ' i+,s, / +Ti+ ' r / comes from symmetry properties of spin and isospin matrix 
elements 

< SfMffaWlSiMi >= (-l) s ^ s f < S f M f \a s m (l)\S t M t >, (A23) 
and the reduced spin and orbital matrix elements are 

< 5/11^(1)11^ > = (-l) 1+9+Si V2 s Si S f | J f | (A24) 

< ly f \\YiK > = ^ZZ^,(C0/0|Z; / ) (A25) 

(A26) 



The radial matrix element 7y ~. ~. is 



7- / \ f°° , 2 [°° , 2 1 L f S f> L 'f S 'f * , \ 

x n^^^J (A27) 
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which can be reduced to a one-dimensional integral over the p-variable 



,2 _ (L\l 



l r°° . \ P ~ [ n' XX 



x xSS A Wxt iisl W ( A28 ) 



r 

where cos a = — and a = a\ = 
ap 



APPENDIX B: 



To calculate the radial formfactors we need the multipole decomposition of the effective 
NN interaction. For this purpose it is convenient to use momentum representation 

V(r pt ,Ppt) = -j—^ J dk exp(-zfc ■ r pt ) V(k,p pt ) (Bl) 

where in V(k,p pt ) the longitudinal v\>{k) = ti T {k) —2t^{k) and transverse tj^(k) = tf T [k) + 
tlj,{k) parts are usually singled out 



V(k,p pt ) = \H^ T (k)a^ ■ of + ^t T LS {k)k xp pt -S- 4{k)S pt {k) >r; ■ r t 

T { S K ) 

= Y.{ t OT^)+M k ){ (T p-k)(^t-k)+t^(k)[(T p X fc] ■ [<T t X fc] 

T 

+ ^t T LS (k)k x p pt ■ 5}r p T • rj (B2) 

Formfactors t l j(k) are Fourier transforms of corresponding forces in coordinate space 

t c ST {k)=An Jo (kr)t c ST (r)r 2 dr, (B3) 
Jo 

roo 

£(k) = 4vr / j 2 (kr)t^(r)r 2 dr, (B4) 
Jo 

poo 

t T LS (k)=Ank hikr^lsirydr. (B5) 
Jo 

With shorthand notations for multipole operators, 

Pisj, m (i) = jiikr^lY^r,) <g> af[ jm (B6) 
P l ij, m (i) = piiknWiri) <g> Vi} jm (B7) 

Ptrni 1 ) = -^(V l pjO j , m { i )) XPi-<7i (B8) 

the multipole decomposition of the NN potential can be written as follows 



JT T 
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2 f°° 

V(r pU ppt) = Y, t p ■ T I ~ I dkk 2 { t 0T (k) (p j0j (p) ■ Pjoj(t)) 

jrp 71 JO 

+ 4(*0 (p|(p) • /#)) + 4(k)( (pf(p) ■ pf(t)) - (p iy (p) • p,i,(t)) ) 

- \t T LS (k) ( (^(p) ■ p j0j ( t )) + (p j0j ( P ) ■ P i f(t)) + (pf( P ) . 4(t)) 



(B9) 



Inserting this decomposition of the iViV-interaction into the expression for the nuclear 
formfactor < J/M/ | St ^pt I JiMi,j a m a >, we obtain |y] formula fll4l) . 

The radial part of the formfactor can be written in the following form: 
a) for excitation of normal parity states, 



l-U (><■>■,, ^-) = E < T t M T TM T \T f M Tf >< T b M Tb TM T \T a M Ta > 



_-\)T b -T a 1 I 

-iU (-||r T |P 

T \ 2 \\ P II 2 , 



x 



.„• s 2 r 00 



J/ vr Jo I 

(la - Jb)jj(kr p ) 



sQOlj 
1 



jjO P ) t 0T (k)p mT (k) - ^LfMPjOi/K*) 



9 ^-(fcrp 
' la 



1 



(BIO) 



+^ [j'O' + 1) - (7b - la) (lb -Ja+ l)]jj(kr p )j^Y 



-^ tT Ls( k )PjljA k ) 



j(j + l )3i(kr P ) 



k 2 r p dr p 



d [(lb - la) (l a + lb + 1) - j(j + 1)] 



2 vO'(j + 1) 



X 



■1, I I''' " 23 
b) for excitation of unnatural parity states 



+ 5 sl 5ijij{kr p ) t T (k)p jljiT (k) 







/•t:U £-) = E < T^TMrlTfM^ >< T b M n TM T \T a M Ta > (^||rj||-) 



xz 



dvp J 

a 9 /"OO 



-2^ + 76 + 1 - j)< 



+S sl Si J+ xj j+1 (kr p ) 
where 



[jtt 

2j + l 



^tmpiAk) - ^i^tmpfAk) 

(Bll) 



4(% r0) - 



t^(k)pf T (k) 



la = (j a m a \L ■ cr\j a m a ) 



2j + l 



I i=l+l 

l ai Ja l a ~ 2 

~(l(i I)) 3a ^a o 



(B12) 



In the formulas above pi s j,T(k) = (JfTf\\ J2t Pisj(t) T t^\\JiTi) is a complex expression con- 
taining spin-angle reduced matrix elements and one-dimensional integrals over radial parts 
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of different components of bound and continuum wave functions and given in Appendix A. 
Other densities are the different linear combinations 



4r(*0 = y ij-r^Pj-iijA k ) + y ^p^P3+u3,A k ) ( B13 ) 

Pf,A k ) = ^ ^p^Pj-Uj,A k ) ~ \j ^J^jPj+llj,T( k ) ( B14 ) 

In radial formfactors we omit the contributions from the current p\jA k ) = 
(JfTf\\Y,tP\j(t) T 't\\JiTi) and spin-current p l j :T {k) = (JfTf\\Y,tP l j(f) T f\\JiTi) densities 
which we did not take into account in calculations. 

The transition densities in coordinate and momentum space are simply connected by 

2 r°° 

pisj,A r ) = - d( i q 2 ji(qr)pisj,T(<i) ( B15 ) 
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